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Abstract 

Flight control of small crew return vehicles during 
atmospheric re en try will be an im p ortan t technology in 
any human space flight missions undertaken in the 
future. The control system p resented in this paper is 
applicable to small crew return vehicles in which 
reaction control system (RCS) t hrus ters are the only 
actuators available for attitude control. The control 
system consists of two modules: (i) the attitude 
contiolicr using the trajectory linearization control (TLC) 
technique, and (ii) the reaction control system (RCS) 
control allocation module usiqg a dynamic table-lookup 
technique. This paper describes the design and 
implementation of the TLC attitude control and the 
dynamic table-lookup RCS control allocation for nominal 
flight along with design verification test results. 

1. INTRODUCTION 

One of the goals of the NASA Orbital Space Plane 
(OSP) program was to provide a small crew transfer 
vehicle for the International Space Station, using only 
technology which had already been matured to a level 
suitable for human flight Two concept vehicles 
proposed to meet this challenge were a small winged 
lifting-body (WLB), and an axial-symmetric capsule 
(ASC) similar to the Apollo reentry module. Both 
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vehicles used only RCS thrusters for. attitude control, 
however die aerodynamic characteristics are drastically 
dif fe rent. The techniques pre se n ted to control these 
vehicles during reentry consist of two modules: (i) the 
attitude controller using the trajectory lin earizati on 
control (TLC) technique, and (ii) die control allocation 
module using lookup table based dynamic inversion and 
dynamic torque matching technique. 

Nonlinear tracking and decoupling control by 
tr aj ectory linearization 1 is one of the advanced 
nonlinear control techniques that can be viewed as the 
ideal gain-scheduling controller designed at every point 
on the flight trajectory. A p rototyp e tr aj ec to ry 
linearization controller (TLC) design for die X-33 
ascent 2 and entry 3 flight phases was developed under 
NASA MSFC Advanced Guidance and Control Program 
and NASA Space Launch Initiative NRA8-30*' 3,6 . For 
each flight phase, the controller has only 36 tunable 
gains that remain constant fur nominal flight, 
dispersions and failure conditions. Initial testing of the 
TLC ascent and entry controllers showed encouraging 
results. During the 2003 NASA MSFC Integration and 
Testing of Advanced Guidance and Control (ITAG&C) 
program 7 , the TLC control algorithm was able to fly the 
most dispersion and failure test cases. In fact, the only 
test cases that the TLC design could not fly were the 
cases that no other algorit h ms could fly at tint time. 

The TLC employs an open-loop controller that 
acquires the nominal guidance trajectory, and a closed- 
loop controller that asymptotically eliminates (stabilizes) 
the guidance command tracking error due to vehicle and 
enviro nmenta l dispersions and vehicle failures. 
Therefore it provides robust stability and performance 
at all stages of flight without interpolation of controller 
gains, thereby eliminating costly controller redesigns 
due to minor airframe alteration or mission 
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reconfiguration. The TLC computes the commanded 
torque symbolically based on the theoretical 
derivations, which am applicable to any feasible 
trajectories and airframes with known mass properties. 
Therefore the design is vehicle scalable and mission 
adaptable. In particular, it should be applicable to the 
WLB and ASC witho ut mo d i alteration. 

To examine and demonstrate the vehicle scalability 
and mission adaptability of the TLC control technique, 
toe TLC control algorithm is applied to both of die 
a fo re m e nti o n e d crew transfer vehicle models in entry 
flight phases with a set of fixed gams for the nominal 
flight conditions. The gam parameters along with initial 
design ve rifica t i on test simulation results are presented 
in this paper. The v e rificati on simulations are performed 
in TDOS, which is a MATLAB/SIMUUNK 6-DOF 
launch vehicle simulation e nvironment , deve lo ped by 
Universal Space Lines LLC under NRA8-30* . 

Another flight controller design challenge for the 
crew return vehicles is the control allocation, which 
realizes foe attitude controller's torque command by 
selecting RCS th r u s te rs, or deflecting a erodyna m ic 
control surfaces if drey are available. The only available 
control effectors for the WLB and ASC concept vehi- 
cles are RCS thrusters with limited fueL The discrete- 
time, discrete-thrust na t u r e of the RCS control not only 
introduces large enure in foe control torque realization, 
but also in troduces high-frequency noises and signifi- 
cant delays in foe control loop. These pro b l ems will 
increase guidance command tracking enure and 
decrease the robustness and stability margin of the 
attitude control. Further compounding the problem is 
foe sensitivity of foe vehicle dynamics and RCS control 
to the offset of vehicle CO. The luge number of 
attainable control moment (AjCM) vectors from all 
combinations of foe RCS thrusters, while providing 
redundancies in case of RCS failures, renders control 
allocation by enumeration of the ACM vectors 
impractical. 

A lookup tabic algorithm is presented to compress 
foe laige amount of ACM vector data into a format 
which can be used to quickly generate a set of RCS 
thruster commands to closely march a commanded 
vehicle body torque. Common practices in control allo- 
cation 9 ' 10 seek to map instantaneous torque commands 
to the vehicle control surfaces and/or RCS thrusters. 
An alternative approach is a dynamic allocation where 
foe torque is allocated by matching a commanded 
torque over a past time interval, or matching a torque 
impulse. In addition, the RCS actuator dynamics will be 
accounted for by * dynamic pseudo-inverse. Preliminary 
control allocation results arc included in this paper. 


Following this introduction, Section 2 pre se nt s an 
overview of the generalized TLC control method 
applicable to reentry space vehicles. The table lookup 
control allocation method is presented in Section 3. 
Preliminary simulation and verification results from 
IDO$ are presented in Section 4. Section 5 concludes 
flic paper with a summary of foe main results and issues 
that need to be addressed by farther research.. 

2. Trajectory linearization Attitude 
Controller 

Consider output tracking by a nonlinear dynamic 
system 

e(f)-/(£(f),*i(t),*(t)) (i) 

•?(*)= 

Where £(<) denotes the State trajectory, ij(t) fo c output 
trajectory, (*{t) the control input, and 0(£) time 
dependent para m eters . Let be the 

nominal state, output trajectories and nominal control 
satisfying 

»?(*)- h^{t),fi(t),0(f)) 

Define the state and output tracking errors, and the 
tracking error control input by 

«(0=€(0-¥(t) 
q(*)-»7(*)-» 10 ) 

The nonlinear tracking error dynamics can then be 
written as 

k(t ) -/(?(*)+£(<), P(t)+Mt), 0(tj) 
-/(!(t), /!(*),*(£)) 

Asymptotic tracking can (hen be achieved by a 2 
Degrccs-of-Freodom (DOF) controller consisting ofi (i) a 
dynamic inverse I/O mapping q ■-» f* of the plant to 
compute foe nominal control jt(f) for any given nominci 
output trajectory q(f), and (ii) a tracking error stabilizing 
control law ft to account for modeling uncertainties, 
disturbances and initial conditions.THe overall 
controller configuration is given in Figure 1 . 
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Figure 1. Ncnlmcw Tracking SySero 
Configuration 


The trajectory linearization tracking controller nacs . 
a linear time-varying (LTV) stabilization controller to 
stabilize the nominal (command) trajectory Jj with a 
linear state feedback control law Ji For the 

WLB and ASC flight cortrollera, the LTV control law 
will be computed using para) lel-difierential (PD)- 
cigen structure assignment A detailed theoretical back- 
ground on nonlinear inversion and PD-dgenstrocture is 
given in 1 - To make this paper self-contained, the main 
design results are su m mariz ed below. 

State variables of the WLB and ASC attitude 
equations of motion in entry and launch short flight are 
the aerodynamic angle vector tf«nd body rate vector ur, 
and flic control input is the body torque vector T, 
which are given, respectively, by 


V 


P 


Tt 

a 

, Um 

q 


Tm 

.0 


r 


T n 


where ip, a, and P are toe bank, angle of attack, and 
sideslip angles; p, q, and r are the roll, pfodb *nd y*w 
rates; and I), T™, T n ■« the roll, pitch and yaw 
moments about the vehicle center of gravity dne to 
airframe and control surface aerodynamics, and RCS 
thrusts. The kinematic equations tj = /„(»*,«) we given 
by. 

y p - [pcos(a) + rsin(o)]/cos(4 + ® 

- t»n(0) [cos$ )cos(y>) - ^aocetl G/ v 
a = q- tan(jS)tpeos(or) -t-rwi(o)] 

+ (cos(d)cos(y) - ZwaxA ]G/V/cos{fl) 
#=p»n((*)-rcos(a) 

+ [cos(d)*m(y>) - Ymxd\Gf v 

where x « the heading angle, d is the flight path angle. 
V is- the magnitude of the relative velocity, G is the 
gravitational constant, Ijaxl and Z t aTC *** 
and normal accelerations in flic Wind Frame, 
respectively. The body rate dynamic equations i* = 


p= 7* g pq+f qr qr + !$(CN p p+Cb r r + 7} > ) (3) 
+ (pN r P + C N r r + Tn ) 
q=J%,f?+ fir r ’ + *pr P r + 9m ( C ^+T m ) 
r - IpqP<l+rqr« r + fl ( C h p + C ^r r + J 0 
+ ( C N,P + C N r r + Tn ) 
where pj*, etc. are coefficients consisting of the 
moments of inertia; 



C U ~ 


9L 

dr 


are 




m 

dq 


„ *N „ 

C AV - Qr 

the dynamic derivatives, and J is the inertia matrix 


J = 



0 

% o 

0 (b 


Owing to the block-triangular form of the equation* 
of tnotion<2) and (3), the plant admits an hmmocta 
loop controller design which greatly simplifies dm PD- 
dgenvaJnc assignment and the nonlinear dynamic 
inversion. However, the validity of this design rests on 
the requirement thsit the timer loop dynamics must be 

sufficiently f»st so as to render it as a singular perturba- 
tion to the outer loop. The trajectory linearization 
controller configuration shown in Figure 1 is applied to 
both the attitude error feedback loop and angular rate 
error feedback loop. 

For attitude angle tracking, the nominal body rate 
«=[p ~q 7] t is given by 

p *= 7pcoi(a) — {*in(y>) — P*xet) *®(a)C?/V (4) 
— i*m(d) cos(5) 
q = S - (cos (ip) - Z^OG/V 
r =* y5sin(5) + [ah(?) ~ Pacccll cos(3)Gf/V 
-Xsm(d)sin(o) 

For body rate tracking, the nominal body torque 
T —\Ti I'm T„] t , assuming Jry •- V 3 °> * 

given by 
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5 \ ~ JxxP+ ( /» — /*|r)ff r (^) 

- 4*(r + qp) - C h v - Cl? 

Tm = fyifl + (/** — IzzVf* 

+ /xrCP 2 -? 2 )~C Mf q 
T n = IxF + (/» - IzxJqP 

+ -Cfi r P~ C N ? 

Integral feedback is employed for disturinmcc accom- 
modati on and robustness to parametric uncertainties. 
The attitude and rate tracking error state variables are 
defined, respectively, by 

I{<P — V > coni)df /(p-Jteot n)dt 

ip — ipixm P—Pctm 

Rot — <*com)cft _ f(.4-<Jam)dt 

T&UR — a — atom *‘^8 q — <kam 

RP-pvxtddL J(r- retm )dt 

0~ fleam t" ~ r com 

where commanded bank angle, fpcom end commanded 
angle of attack, acorn are given by the eloserMoop entry 
guidance, and the sideslip angle — 0° must be 
tigh tly regulated. The proportional-integral (PI) 
feedback control law for the attitude loop tracking error 
is given by «i — -Xi(t)tfc^, where the gain matrix 
K\ (t) is computed symbolically as 

*l(t) = («) 

*H1 *112 0 ^ *116 *116 

0 *122 ®12l C*I22 0 *126 

*131 *132 0 ~P *136 *136 

where 

*111 = cos(a)oiii 

*112 = cos(a)aii2 + sin(a)cos(p)G/V 
*115 = san(a)ai31 
*116 — sin (5) ai32 

- (cos(p) - ^wseci) cos(a)G/V 
*122 «= -sin &)G/V 

*126 =— [pco*(5) + r sin (a)] 

*231= sin (a) an i 

*132= sin(o)aii2 -cos(5)cos(p)G/V 
*135 = — cos(a)r*j3i 
fejje = -cos(S)ai32 

— [cos (ip) — i?aocell (a)G/ V 

Similarly, tile Pi feedback control law for the rats loop 
tracking error is given bytij *= -K 2 (t)v ing- The gam 
matrix K% (t) is computed symbolically by 


K 2 (t)~ , (7) 

*211 *212 0 *214 *216 *216 

0 *222 *223 *224 0 *220 

*131 *132 0 *234 *236 *236 

where 

*211 = I x* arm 

*2)2 =/*e(i^? + «Wll) - /»'^5 + C!£ > 

*214 *= /ss(/pjP + /flr r ) “ /p-f) 

*215 =>—/** 0331 

*216 *= /**/£?- 1x^7+0212)+^ 

*222 — + /pr r ) 

*223 =/w q 221 
*224 = Cfd n + /w0222 

*226 = i*|r(2 tr r + fj^p) 

*231 - -/x*a31) 

*332 = IsftpqQ Zxr(Z pq 7 **212) + Cff f 

*234 = Tu(I^P+ fi/r r ) ~ V + tffr r ) 

*235 = Ixz »231 

*236 - Ia{IjrV + 0232) - !**/£■? + C N r 

In the above gain matrices K\(t) and Ki{t), tire 
parameters Oi//(t), *=1,2, y=l,2,3, 1= 1,2 are 
obtained from tire closed-loop quadratic PD-dgenvahres 

mo— ( a# * V 1 (8) 

for Qj ^ 1 , with constant damping Qj and time-varying 
bandwidth uhij(t) as follows 

«yi (<)*“&/*) & 

^ nij (f ) 

aij2(t) = 2(ijutoij(t)- 

Closed-loop stability with respect to the equations of 
motion is established far > 0 and Wnij(f )>0, for 
all t > 0. As a eonstraint of tire singular perturbation 
method, the inner loop bandwidth u>n2 j(t) must be 
sufficiently greater than the outer-loop bandwidth 
to ensure stability. In addition, fa =1.0 
should be avoided. The time-varying bandwidth wn)y (t) 
can be used for adaptation to vehicle and environmental 
dispersions and control effector failures. 

The reference signal for tire outer loop is the 
commanded attitude fjtom = 7- The inner loop reference 
is the commanded body rate u> C om = ui + «i . The 
controller output is the commanded torque 
2’ com = T+ 1*2» which is the input to the control 
allocation unit. The commanded torque is realized via 
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die neural network based control allocation unit 
discussed in the following section. 

Due to the integral (high-gain) control, loss of 
control authority due to excessive disturbances such as 
wind gust, vehicle dispersions such as overly estimated 
control eff ect iv eness, and control effector failure such 
as loss of one or more RCS t hr uster s , will result in 
torque command saturation. This in e ffe ct increases the 
level of singular p er turba tions to the outer control 
loops, which may result in performance degradation or 
even loss of stability. By the singular perturbation 
theory, reducing die outer loop bandwidth accordingly 
can ensure stability and p er fo r ma nce, although the 
overall response time will be longer. The time-varying 
bandwidth capability of the PD-eigenvahje assignment 
technique allows the vehicle to adapt to severe 
dispersions and partial control effector futures without 
explicit fault identification. This is known as direct fault 
tolerant control, and has been demonstrated to be very 
effective under the NRAK-30 Program 1 1 . 

The overall flight controller consists of fi PI 
controllers, one for each channel hi the inner and outer 
loops (note that the outer loop is effectively controlled 
by a proportional-integral-derivative (PID) controller, as 
the inner loop serves as a deri va tive feedback). Each 
channel Has a nominal controller c ont aini ng a pseudo- 
differentiator. With a constant rate time- varying band- 
width each channel in a loop has 6 constant 

design parameters to be tuned for stability and 
performance. The total number of tunable design 
parameters is 36, which include: Wnit.difT, fc= 1 .2,3, 
which arc the bandwidth of die outer loop pseudo- 
differentiator for the roll, pitch and yaw channel, respec- 
tively; ft*, which is the damping ratio and u^i^rnh,, 
Hilt, max which arc foe lower and upper limits of foe 
closed-loop bandwidth for tire attitude error dynamics; 
“blMoc «id Whjkfac which are tire decrement and 
increment of the outer loop bandwidth per control cycle 
when any of the control moment command exceeds or 
recedes the m ax i mum attainable moment in that channel 
(torque command saturation), respectively. The inner 
loop counterparts C 2k. WnJMiin. ‘‘b 2k, max, 

s »id i*,j*,inc, k *■1,2,3, are defined similarly. 
These constant design parameters, once tuned, arc used 
for ah nominal and abort flight conditions without 
alteration. 

3. Control Allocation Design 

The commanded torque T^m generated by the 
trajectory linearization attitude controller is realized via a 
control allocation unit that produces the corresponding 
RCS firing command 6, which is an AT-dimcnsional, 


binary-valued vector, where N is foe number of 
thrusters, the ih entry Si corresponds to the tth ' 
thruster 4 * on-off command with 6* »* 1 for t hrus t e r on, 
and. 6i = 0 thruster off. It is a common practice to 
formulate foe control allocation problem as finding a S 
such that 

T = B6 *» Tco* (10) 

where B is foe thrust-to-moment table for all RCS 
t hru ste r s, or the control effectiveness matri x, which can 
be obtained from a given vehicle RCS co n fig u ration. 
The control allocation, or the RCS thruster selection 
logic design thus amounts to finding an inverse 
mapping 3^ -* 6. 

The WLB and the ASC both employ only RCS 
control effectors. For a vehicle configuration with N 
RCS thrusters, the total number of attainable control 
moments is 2^, which constitutes the attainable 
moment set in the roll-pitch-yaw moment space. To 
ensure effective attitude control and guidance co mman d 
tracking, and to mi n i mize RCS fuel consumption, it is 
desirable that the attainable momenta be distrib u ted 
evenly in the moment space. 

The range of moments attainable from the 16 RCS 
thrusters on the WLB are shown in figures 2-4, the 
moment Tange attainable from the 12 thrusters an foe 
ASC are shown in figures 5-7. These figures show the 
entire set of attainable moments in each channel from 
least to greatest magnitude. The segments with gaps 
betw e en adjacent points on the lines indicate 
unattainable moments within the range. The rotl 
moment sets for both configurations indicate a number 
of roll moments that are unattainable by any 
combination of thruster firings. The ASC has 
particularly uneven roil moment distribution, which 
greatly reduces foe controllability of the vehicle. 

The distribution of moments in one channel with 
respect to another channel is important in determining 
foe ability of a control allocation algorithm to meet a 
large combination of eommandod roll, pitch, and yaw 
moments. The distribution of moments in the moment 
space are shown in figures 8-10 the WLB, and figures 
11-13 for the ASC. A dense, even distribution such as 
tire pitch vs. yaw moment plot for the WLB in figure 10 
is desirable because it indicates that almost any 
combination of pitch and yaw momenta can be realized. 
On the other hand, a distribution such as the roll vs. 
yaw plot in figure 12 for foe ASC indicates that lor any 
desired yaw moment, there are few commending roll 
moments, which also reduces foe controllability of foe 
vehicle about foe roll axis. 

Jt is noted that at high angles of attack, both 
sideslip regulation and bank command tracking require 
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coordinated roll-yaw control. Tims the attainable 
stability margin and guidance command tracking 
performance can be improved significantly if the 
orientation and location of the thrusters can be adjusted 
within tire vehicle configuration constraints so as to 
improve die distribution of the attainable roll-yaw 
moments. In fact, the d istrib u tion of attainable moments 
from a thruster c o n fig ur a ti on can be significantly 
im pr ov ed by making flight adjustments in die plume 
directions of a few t hrus ters. Given die impact of 
moment di s trib ut i o n s on the controllability of a vehicle, 
it would be prudent to address this issue early in the 
vehicle design stage. 

Due to the binary-valued RCS control e ffe ctor 
vector 6, the inverse mapping cannot be constructed 
using a pseudo-inverse of the B matrix in equation (10). 
Thruster allocation by table lookup is a simple algorithm 
capable of finding a set of RCS thruster commands that 
produces a body torque which closely matches the 
commanded torque. The lookup table method is die 
most accurate way to allocate control to the RCS 
thrusters daring nominal flight. Since the majority of the 
computations ate done off-line, the lookup table can be 
configured to trade off accuracy for other 
considerations such a comp u t a ti o nal power, memory 
limitations, fuel consumption, etc. The drawback is that 
in order to accommodate off-oomirwl event* such as 
thruster failures, a di ffer ent lookup table would be 
needed for each event, which may not be feasible with 
limited onboard com putati onal resources. 

The lookup table is constr uc ted by first discretize 
the range of attainable control moments along the toll, 
pitch, and yaw axes with a denied resolution. The 
discretized moment values along cadi axis are ordered 
monotonically increasing to produce a roll, pitch, and 
yaw moment analog-to-dighal conversion (ADC) vector, 
respectively, so that each discrete interval along the roll, 
pitch, and yaw moment axes can be indexed with an 
integer, and each discrete (rectangular) cell in the 
attainable moment set can be indexed with a 3- 
dimensional index vector that serves as die indices for 
the lookup table. 

Each attainable moment is represented by an N- 
digh binary number co rresp o nding to the RCS thruster 
firing combination 6. The lookup tabic is then construct- 
ed by choosing one attainable moment from each cell. 
Each cell of the table is thus filled with an integer 
corresponding to die binary number. If there are more 
than one attainable moment in a cell, the one that uses 
the least number of thrusters will be chosen. If a cell is 
empty, then an attainable moment from a nearby cell that 
is in a certain sense closest will be chosen. 

A trade off between accuracy and computational 
constraints can be met by adjusting the resolution at 


which the moment ranges are discretized. Since the 
discretized ADC vectors need not be the same length, a 
certain channel may be favored by discretizing the 
moment set of that channel at a higher resolution. The 
method of determining the closest cell to an empty cell 
can also be tailored to favor a certain channel by only 
searching along the other two axes. A logic algorithm 
could p ro duc e one of several tables with different 
perf orm a nce characteristics based on ament vehicle 
conditions. 

Lookup tables were created for the WUB and ASC 
thruster configurations. The tables for both used 60 
partitions along each axis which resulted in lookup 
tables with 60 3 elements. The lookup tables were then 
used to find thruster commands to meet the torque 
commands from a typical entry torque command profile. 
The body torque induced by the thruster commands 
found from the lookup tables are plotted over the 
commanded torque in figures 14-16 for the WLB, and 
figures 17-19 for the ASC. 

4. Simulation Test Results 

At this time, preliminary tost results for the TLC 
controller are shown in Figures 20-31. The test 
simulations were p erform ed in TDOS, which is a 
MATLAB/SIMULINK 6-DOF bunch vehicle simulation 
e n v ir o nm ent. Figures 20-22 and 27-28 show e n ter -loop 
tracking performance for flic WLB and ASC, 

respectively. Figures 26-28 and 29-31 show timer-loop 
tracking perfor man ce for the WLB and ASC, 

respectively. Although the two vehicle models differ 
drastically, in these initial tests the TLC controller* for 
these two vehicles used the same gains- While further 
tuning with increased fidelity of the models will likely 
lead to different guns, the current results have 
demonstrated the vehicle scalability and mission 
adaptability of the TLC control technique. 

5. Summary and Conclnsions 

In this paper we have demonstrated the vehicle 
scalability and mission adaptability of die trajectory 
linearization control (TLC) technique by applying it to 
the design of entry flight controllers for two drastically 
different small reentry vehicle— a winged lifting body 
configuration and an axial-symmetric capsule con- 
figuration. Initial test results showed satisfactory 
performance with the same constant controller gains. 
Moreover, each flight controller only uses 18 constant 
controller gains when no adaptation is employed, and 
with 36 constant controller gains direct fault tolerance 
adaptation can be implemented. 
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In addition to the TLC attitude controller, a table- 
look up RCS control allocation design method is 
presented. The method is simple yet effective, as has 
been demonstrated with preliminary test remits. The 
design presented herein is only the first step of a 
dynamic control allocation technique. Additional week 
is being p er formed to implement the dynamic table- 
lookup method which provides best match of the 
commanded torque over a recent past interval of time, 
rather than the matching the instantaneous torque 
command. 
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